Reconstruction of Rb-Rb inter-atomic potential from ultracold Bose-gas collision 
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Scattering phase shifts obtained from ^'^Rb Bose-gas collision experiments are used to reconstruct 
effective potentials resulting, self-consistently, in the same scattering events observed in the ex- 
periments at a particular energy. We have found that the interaction strength close to the origin 
suddenly changes from repulsion to attraction when the collision energy crosses, from below, the 
I = 2 shape resonance position nt E ~ 275 fiK. This observation may be utilized in outlining future 
Bose-gas collision experiments. 
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Introduction Scattering and collision processes have 
always been one of the most important tools providing 
information on interaction between quantum systems or 
particles at all lengths and energy scales. In the quantum 
regime, we usually cannot perform direct experiments 
with atoms, electrons or other particles. What we can 
do is merely to reconstruct scattering events and thereby 
the interaction properties of the constituents from exper- 
imentally observed quantities of the collision. At high 
energy, particle-colliders serve to probe Nature at the 
smallest distances. At the lower end of the energy scale, 
we can cool down atoms to the ultracold regime (< 1/iK) 
P, 0| and form a Bose-Einstein condensate (BEG) ^ to 
explore the low-energy properties of atomic interaction. 

At low temperatures, collisions play a pivotal role, 
mainly in affecting the static and dynamic properties 
of the condensate: stability, lifetime and thermalization 
rate. Furthermore, it is much simpler to describe the 
dominant binary collisions theoretically than at higher 
energies, because inelastic processes are usually negligi- 
ble. In early BEG experiments, the s partial wave and the 
associated scattering length Og [H have been sufficient to 
characterize the low-energy properties of the condensate. 
However, to produce molecular condensate from alkali 
atoms, or to tune the interaction strength via Feshbach 
resonances, or to discover effects beyond the mean- field 
theory, we need to know the interatomic potential be- 
tween atoms more precisely. Both theoretical calcula- 
tions [1, 0, 01 and experimental efforts 113, Ei have 



been devoted to this aim, e.g., the authors of Ref. [ll| 
presented high-precision fluorescence data and accurately 
fitted parameters of the electronic ground state of Rb2. 
Theoretical calculations usually involve some experimen- 
tal quantities, like scattering lengths (6|, positions of Fes- 
hbach resonance 0, or Raman transition rates [^. 



These results are then compared with other data, and if 
the correspondence is not satisfactory, the potential used 
in the calculation is iteratively adjusted. 

Our approach is conceptually different. We want to 
characterize or even reconstruct the inter-atomic poten- 
tial from scattering phase shifts, rji, by using the in- 
verse scattering method. Inverse scattering theories (see 
I3I . 14 1 and references therein) provide two types of po- 
tentials. One type is derived from phase shifts measured 
at a particular angular momentum but for all energies. 
The other type of potentials is calculated from all partial 
wave phase shifts corresponding to a fixed energy. The 
underlying idea of both types of inverse scattering theory 
is the assumption that there exists an "effective" spheri- 
cal potential which is the cause of the observed scattering 
events. In this Letter, we employ the fixed-energy inverse 
scattering theory and use experimental phase shifts rj^^^ 
derived from Bose-gas collision experiments p,5i, ld\. Be- 
cause the experimental phase shifts are largely riddled 
with errors, we rely upon the phase shifts rjf 77™^ of 
Julienne's coupled channel calculation [20J, by which the 
BEG collision experiment [l3| has been interpreted. No- 
tice that these phase shifts stem from a coupled-channel 
calculation and reproduce the well-known d-resonance in 
the I — 2 partial wave at collision energy E w 275 fiK. Al- 
though this resonance is considered as a shape resonance 
of the triplet ^^Rb - ^^Rb potential [2l[, without its 
inclusion one cannot explain [8| important coupled chan- 
nels effects like tunability of scattering length, which is a 
typical Feshbach resonance effect. We demonstrate that 
the coUisional phase shifts derived from experiments also 
contain such information, if we consider them at fixed 
energy but all partial waves. Note that this considera- 
tion is precisely the situation which occurs in collision 
experiments. 
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In the following, we briefly review the Bose-gas colli- 
sion experiments |15l.[l6j which provide the input data for 
the inversion. Thereafter we collect all the the formulas 
necessary to reconstruct potentials from phase shift data. 
In discussing the results, first we present the ^''Rb-^^Rb 
inversion potentials extracted from measurements carried 
out in the low collision energy domain 100-200 /iK, then 
we exhibit these effective interactions in the resonance 
region around 275 ^K. Here a sudden change from repul- 
sion to attraction of potentials can be observed and in- 
terpreted as the coupled channel effect mentioned above. 
Finally, we show that well beyond the resonance region, 
between 600-1200 where only calculated phase shifts 
are available, the potentials appear to be independent of 
energy and retain their attractive character. 

Review of BEC collision experiments In the exper- 
iment of Ref. [ll], *^Rb atoms were prepared in the 
\F ~ 2,mp ~ 2) hyperfine state and precooled to about 
12 /iK. Afterward the cloud of atoms was adiabatically 
split into two pieces by an emerging potential barrier in 
the middle of the trapping potential. Thereby the atom 
cloud experienced a double well potential and was di- 
vided into two parts separated by 4.3 mm. The barrier 
height was several times higher than the chemical po- 
tential, therefore the clouds were isolated and could be 
further cooled down to 225 nK. After such a prepara- 
tion, the double-well potential was ramped down to a 
single-well harmonic potential characterized by the an- 
gular frequencies = 27r x 155 Hz and Wz = 27r x 12 
Hz, and was kept constant during the measurement. As 
a result, the stacks of atoms were accelerating towards 
each other and finally they collided with a relative ve- 
locity V. The collision energy range i?kin/^B = ^v'^/2kB 
was between 87 and 553 /iK calculated in the centre-of- 
mass frame. Due to symmetrical collision of identical 
bosons, the p partial wave was prohibited. Following the 
impact, the atoms were scattered and were moving in the 
trap until their maximum extension. At this moment, an 
absorption image was taken with a resonant light shone 
onto the clouds, perpendicular to the scattering axis. Al- 
though this image was a two-dimensional projection of a 
three-dimensional density, it was possible to reconstruct 
[ist the full tomographical information taking into ac- 
count the cylindrical symmetry. 

In a similar experiment [l^l , it has been demonstrated 
that it is not necessary to keep the atoms in the trap 
all the time after the collision. The analysis of the ob- 
served scattering halo provides all the data. Calculating 
the density in small angular sectors yields us the angular 
scattering distribution of the halo, which is directly pro- 
portional to the differential cross section a{9). Taking 
into account only the 1=0, 2 partial waves and fitting the 
analytical expression 
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to the measured angular distribution finally provides the 
experimental phase shifts ?7q'''^ and rf^"^. These latter 
phase shifts constitute the input for our inverse calcula- 
tions. 

Inverse scattering formulas We employ the fixed en- 
ergy inverse scattering method of Cox and Thompson 
(CT) [l3| in order to derive model independent potentials 
from the given phase shifts 3, corresponding to the 
measurements. The CT inversion method has a number 
of useful properties 3, we may work with a finite 
set of N experimental phase shifts and obtain inversion 
potentials of non-zero first momentum, J rV{r)dr ^ 0, 
and of finite value at the origin. Let us denote the set of 
physical angular momenta I by S and the set of unknown 
'shifted' angular momenta L by T. The latter has to be 
determined from the phase shifts. T contains the same 
number TV of elements as S, and the sets 5" and T are 
disjoint. The CT method leads to a system of non-linear 
equations 
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(1) 



in which the input scattering phase shifts {5i) determine 
the 'shifted' reactance matrix elements defined as 

^ ML{M-')^,e^^^'-'^-/', (2) 

with the square matrices 

N'\ _ 1 Jsin((? -i)7r/2) 

M 



LiL + 1) - l{l + l)\cos {{I -L)7r/2) J, ^ 



(3) 



containing the unknown L— values. Once the set T is 
determined by solving the highly nonlinear equation ([1]) 
for the Ls, we calculate coefficient functions Ai^{r) using 
the system of linear equations 

l{l + l)-L{L + l) -"'^"^' 

where ji^ and mean the spherical Bessel and Neumann 
functions, respectively, and VF[a, 6] = ab' ~ a'h denotes 
the Wronski determinant. Next, we compute the trans- 
formation kernel 
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from which the inversion potential is obtained as 

2 d K{r,r) 



V{r) 



■ dr 



(5) 



(6) 



Results and discussion In Fig. [T] we present the re- 
sults for the energy range between 100-200 /iK, which 
lies below the characteristic I = 2 resonance position of 
the *^Rb - *^Rb triplet scattering. The inverse poten- 
tials are repulsive at smaller distances and oscillatory for 
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larger relative coordinates. The strength of the repulsion 
is approximately proportional to the scattering energy, 
as is the attractive first minimum, the position of which 
gets smaller values from 25 nm to 18 nm. 

Not unexpectedly, the inversion potentials reproduce 
the input phase shifts well within the considered energy 
region. This fact is demonstrated in Fig. [21 where both 
input and output phase shift values are exhibited. The 
phase shift reproduction is rather accurate for the partial 
waves I — 0, 2, 4: involved in the inversion procedure. But 
it must be so because the only control over the potentials 
is the reproduction of the initial data, since there is no 
free parameter in the inversion calculation which is, in 
principle, unique Numerical uncertainty sometimes 
may lead to false potentials which can be recognized on 
false reproduction of the input phase shifts. 

It is therefore important that the excellent reproduc- 
tion of input phase shifts proceeds further on up to 400 
/iK, that is throughout the whole I = 2 resonance region. 
The corresponding inversion potentials are exhibited in 
Fig. [31 and we observe that these 'on-resonance' po- 
tentials produce an abrupt change of potential strength 
V{0) from repulsion to attraction as the collision energy 
crosses the I — 2 resonance position at « 275 /iK. This 
behavior is quite unexpected in view of the input phase 
shifts depicted by lines in Fig.[2l which exhibits a smooth 
behavior in the s— wave and a smooth shape resonance 
in the d— wave. Because of this quite unusual behavior of 
the inversion potential we have performed another con- 
trol test of the results besides the comparison of input 
and output phase shifts mentioned above. We have in- 
verted the original phase shifts of Ref. at two collision 
energies E = 203 and 447 /iK. The values of these mea- 
sured phase shifts differ slightly from the ones obtained 
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FIG. 1: Inverse potentials as functions of interatomic distance 
obtained from the Z = 0, 2, 4 scattering phase shifts within 
the range E= 102-203 /iK, below the d-resonance. Inset shows 
the central amplitude of the inverse potential as a function of 
collisional energy. 



by the coupled channel calculation 15|. Therefore we 
can test at the same time both the stability of the inver- 
sion procedure and the sensitivity of the potential against 
small errors in phase shifts. The inversion potentials ob- 
tained from the data of Ref. [IB] also resulted in the sud- 
den change and were practically the same as the ones 
calculated from phase shifts of Ref. . 

Let us look at the results from another viewpoint too. 
During the collision process the energy is fixed by the ex- 
perimental device and the colliding atoms are approach- 
ing each other from asymptotic distances to the smallest 
ranges. This is described in terms of partial waves. How 
many partial waves are involved depends on the details 
of the interaction of the colliding atoms, just as the very 
numerical values of the partial wave phase shifts do de- 
pend on that. In these Bose-gas experiments, it has been 
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FIG. 2: Phase shifts corresponding to the first three allowed 
partial waves {I = 0, 2, 4) as functions of the collisional energy. 
Solid lines represent the original input data. Symbols x, ▲ 
and □ stand for the phase shifts calculated from the inverse 
potentials shown in Figs. [1] [3] and |4l 
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FIG. 3: Inverse potentials as functions of interatomic distance 
obtained from the i = 0, 2, 4 scattering phase shifts around 
the d-resonance, J5=200-400 /iK. 
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found that the first three {I = 0, 2, 4) phase shifts have 
utmost importance at each fixed energy and these phase 
shifts contain all the information about the interaction of 
the colliding *^Rb gas particles. We have put these three 
phase shifts belonging to each particular energy into the 
inversion procedure, which provided potentials that pro- 
duced the same scattering events as those observed in 
the experiments. Since the inversion procedure is unique 
[13] , we may assume that the inversion potentials are the 
effective potentials which govern the collisions. 

In Fig. m we show the high-energy (well above the res- 
onance) inversion potentials which do not depend too 
much on the energy but retain their attractive character. 



Since the input data stem entirely from experimen- 
tally confirmed data and the inverse procedure provides 
unique results, we do expect that the sudden change of 
strength of the potentials can be observed and utilized 
in future Bose-gas experiments. Moreover, the method 
of inverse scattering may gain ground in the analysis of 
low-energy collisions as it does in high energy physics. 

This work was supported by the Hungarian Scientific 
Research Fund, under contracts OTKA-T47035, T49571 
and the MTA-DFG grant (436 UNG 113/158). D. S. ac- 
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of the Royal Society of New Zealand. 
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FIG. 4: Inverse potentials as functions of interatomic distance 
obtained from the / = 0, 2, 4 scattering phase shifts within the 
range -5=606-951 ^K, above the d-resonance. Inset depicts 
the moduli of central amplitude of the inverse potentials as a 
function of collisional energy. 

Summary The inversion potentials shown in Figs. [U 
|3] and |4] can be viewed as 'effective' interactions which 
characterize the *^Rb interaction in the considered en- 
ergy range, although these Bose-gas collision experiments 
can also be described by sophisticated coupled channels 
calculation where the 'potential' is represented by a po- 
tential matrix. One or more elements of these potential 
matrix may become dominant over the others at different 
energies. In case of the d-resonance, the triplet potential 
plays a governing role in explaining the observed cross 
section. In other cases, channel couplings may destroy 
the dominant effects of a particular element of the po- 
tential matrix. This might be the reason of the change 
in nature of our potentials as we vary the collisional en- 
ergy, since even a weak interaction of different elements 
of the potential matrix can cause resonances (e.g., shape- 
or Feshbach-resonance) . Our effective potentials thus en- 
able intuitive interpretation of the interaction between 
colliding ^^Rb atoms and can be considered as a one- 
channel mapping or a local and energy-dependent equiv- 
alent form of the whole interaction matrix. 
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